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We report on a multiband microscopic theory of many-exciton complexes in self-assembled quan- 
tum dots. The single particle states are obtained by three methods: single-band effective-mass 
approximation, the multiband k -p method, and the tight-binding method. The electronic structure 
calculations are coupled with strain calculations via Bir-Pikus Hamiltonian. The many-body wave 
functions of N electrons and N valence holes are expanded in the basis of Slater determinants. The 
Coulomb matrix elements are evaluated using statically screened interaction for the three different 
sets of single particle states and the correlated iV-exciton states are obtained by the configuration 
interaction method. The theory is applied to the excitonic recombination spectrum in InAs/GaAs 
self-assembled quantum dots. The results of the single-band effective-mass approximation are suc- 
cessfully compared with those obtained by using the of k ■ p and tight-binding methods. 
PACS numbers: 71.35.Cc, 73.21.La, 78.67.Hc, 71.15.-m 



I. INTRODUCTION 

Semiconductor self-assembled quantum dots 
(SAQDs)0, S i are islands of one semiconductor, 
e.g. InAs, in a host matrix of another semiconductor, 
e.g. GaAs. The elementary excitations, electrons and 
holes, are believed to be confined in all three dimensions 
by the band gap difference between island and matrix 
materials. The picture of electrons and holes as confined 
elementary excitations with effective mass, interacting 
via Coulomb interactions has been successfully applied 
toward the explanation of many experimental results 
SUBHEllHini]- It is important to establish 
to what extent the effective-mass picture is applicable 
to the description of electronic states of self-assembled 
quantum dots by a systematic comparison of different 
approaches. The self-assembled quantum dots plus the 
surrounding barrier material contain millions of atoms 
and the density functional oh initio calculations are 
not possible yet. Hence in this work we compare two 
simplified approaches, the multiband k ■ p method and 
the tight-binding method with the predictions of the 
effective-mass calculations. The multiband k ■ p method 
[H [H Q E El [H E3 accounts for the proper 
structure of the valence band, including heavy, light and 
spin split-off hole bands. It is however limited to the top 
of the valence band, does not account for the atomistic 
character of the interfaces between the dot and barrier 
material, and is expected to break down as the size of 
the nanostructure decreases. The atomistic structure of 
the nanostructure is captured in either the tight-binding 
[l9l l2fj| | or pseudopotential approaches as developed by 
Zunger and co-workers |2l|. The tight-binding approach 
chosen here is the effective bond orbital model (EBOM) 
[22I I23I |24|. a version of sp 3 tight-binding models. The 
advantage of EBOM is that it extrapolates to the k ■ p 
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approach making a direct comparison possible. The 
disadvantage of EBOM is that it misses the lack of 
inversion symmetry of zincblende structures. 

The single particle energy levels are not measured di- 
rectly. What is measured in, e.g. optical experiments, 
is the emission from self-assembled quantum dots as a 
function of the excitation power, or the number of elec- 
trons and holes in the dot. The electrons and holes 
interact and form multi-exciton complexes. Emission 
from multi-exciton complexes has been measured by a 
number of groups ElllHEESESlHil- 
The higher the pumping intensity is, the more excitons 
are involved, thus the emission from higher excited elec- 
tron and hole states can be observed. The multi-exciton 
emission spectra have been interpreted using quantum 
mechanical methods such as the Hartree-Fock method 
and the configuration interaction method (CI) (U, in 
which the multi-exciton complex states are constructed 
from single-particle states of the system. It is a chal- 
lenge to combine realistic single-particle states calcu- 
lated for the million-atom structures with these quan- 
tum mechanical methods. A number of theoretical ap- 
proaches have been proposed to address this issue, such 
as combining multiband k • p single- par ticle states with 
self-consistent Hartree-Fock method Hal and combining 
single-band effective-mass 0, l33t 134 ^3 or micro- 
scopic pseudopotential wave functions |36l |3jj . 

In this paper, we use a general approach which com- 
bines different multiband calculations of single-particle 
states with the CI method for the calculation of multi- 
exciton states. By using single-particle states obtained 
from the single-band effective-mass approximation, the 
multiband k • p method and the atomistic tight-binding- 
like method, we are able to compare the multi-exciton 
emission spectra obtained from different single-particle 
states and determine both the validity of the effective- 
mass approximation as well as the validity of multi- 
exciton emission spectra as fingerprints of electronic 
structure of quantum dots. 
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II. SINGLE-PARTICLE CALCULATION 

The single particle calculations for self-assembled 
quantum dots started with the effective-mass calculations 
which related shape and size of the dots to the single 
particle energy levels |3S| . As experimental information 
accumulated, more sophisticated approaches were devel- 
oped, such as single-band effective-mass method cou- 
pled with strain calculation (3, 0, , eight-band k ■ p 
method 13 13 E3 > tight-binding methods [3 E| 
and the empirical pseudopotential method |2l| . In the 
following, we briefly describe the single-band effective- 
mass method, the eight-band k • p method and EBOM 
for the calculation of single-particle states in SAQDs. 




A. Effective-mass single-particle states 

Here we use a single-band model with anisotropic effec- 
tive masses of electrons and holes treated as adjustable 
parameters. The Hamiltonians, including strain, read 
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where H s = e xx + e yy + e zz and B s = e zz - \(e xx + e yy ) 
is the hydrostatic and biaxial strain component, respec- 
tively, V§ and V b h a are the potentials from the band offsets 
between island (InAs) and matrix (GaAs) material. a c , 
a v , and b are the deformation potential parameters that 
are also used in the multiband k • p method and EBOM. 



B. Eight-band k ■ p single-particle states 

The eight-band k • p method uses eight Bloch functions 
at the r point of the Brillouin zone as basis functions to 
describe electron states with finite wavevector. As the 
lateral size of SAQDs is usually much larger than the 
lattice constant, it has been widely used in the calculation 
of confined electron states in SAQDs 0, Q, 0. In 
general, the multiband k • p Hamiltonian can be written 
as 
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where E& is the matrix for the band offsets, A's, B's, 
and C's are the coefficient matrices [l2j]. By using the 
deformation potential theory, an additional part H s |12| , 
which has a similar form as -ffk p , can be added to take 
into account the effects of the strain. 



FIG. 1: Band structure of GaAs described by the tight- 
binding-like effective bond-orbital method (solid lines) and 
the eight-band k • p method (dash lines) . 



Fig. ^plots the energy dispersion of GaAs bands (bro- 
ken lines) calculated by the eight-band k-p method. Note 
a spurious crossing between valence bands at wavevectors 
k halfway between the T and X points. This crossing 
may result in spurious valence band states in SAQDs. 
The problem can be artificially removed by adding addi- 
tional terms prop ortional to k into the eight-band k • p 
Hamiltonian [171 ]. 



C. Tight-binding single-particle states 

By using the same number of basis functions as the 
eight-band k ■ p method, EBOM is a sp 3 tight-binding 
method based on an effective fee lattice |3|, i.e., a pair 
of cation and anion in a zinc-blende lattice is treated as 
a single super-atom. The Hamiltonian is given by 

CRa\H EB \'R / a / ) = E p 6 RR 'S aa ' + <5r-r/,t- 

T 

{E xy T a T a ,(l - 5 aa >) + [E xx t 2 + E zz (l - Tl)]6 aa >}, 
(Rs\H EB \R's') = E S S RR , + E ss "R-R'.n 



(R,s\Heb\R-'o>) = y^^E sp T a 8 R - R ,, T , 



(3) 



where |Ra) denotes an orbital a located at site R. E s , 
E p , E ss , E sp , E xy , E xx , and E zz are parameters that are 
chosen to reproduce the conduction-band effective mass, 
band gap, spin-split energy, and Luttinger parameters. 

Fig. ^ plots the energy dispersion of GaAs bands 
(solid lines) calculated by EBOM. As EBOM uses the 
same Luttinger parameters as the k • p theory does, 
both approaches give the same band structure near the 
r point. However, the spurious crossing from the eight- 
band k-p theory is not found in the EBOM's band struc- 
ture. Here, we adopted the two-center approximation in 
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the parametrization |2^| instead of fitting the energy sep- 
aration between the heavy-hole and the light-hole band 
at X point 22]. 

Compared with the single-band calculation, the multi- 
band methods gives more realistic confined states in 
SAQDs. These states lack the symmetries, such as an- 
gular momentum and spin, which are usually preserved 
in the single-band calculation. These symmetries are im- 
portant in the CI calculation because they substantially 
reduce the total number of configurations. In general, 
the multiband eigenstates of SAQDs do not have any 
spatial symmetry due to the effects of shear strain even 
for a disklike or lens-shape dots that have circular sym- 
metry. In addition, the total spin S and its projection 
S z are not conserved due to the spin-orbit interaction. 
However, we will show that the multiband single-particle 
states are polarized and the polarization can be used to 
define quasi-spin. 



III. FORMULATION OF THE MULTI-EXCITON 
PROBLEM 

Our structure contains millions of atoms and hence 
millions of electrons. As long as the total system contains 
an energy gap and a well defined ground state wavefunc- 
tion, the intractable million electron problem can be re- 
placed by a much smaller problem of pairs of excitations 
in the form of quasi-electrons and quasi- holes. Formally, 
any electron state can be expanded in terms of increasing 
number of pairs of excitations, 

* = * + £>r*r + E c T^T + • ••> ( 4 ) 

where is the Hartree-Fock ground state with all va- 
lence states occupied and conduction band empty. '5™ 
is an excited state formed by removing an electron from 
the state i in the valence band and creating a "hole", 
and moving it to the state m in the conduction band, 
creating an "electron". <3>™ is a doubly excited state 
containing two "electrons" and two "holes", and so on. 
The number of electron-hole pairs is in principle not con- 
served and this expansion can be used to describe all ex- 
cited states However, in semiconductors, the differ- 
ence between the kinetic energies of different numbers of 
pair excitations is proportional to the band gap, which is 
much larger than the Coulomb interaction mixing them. 
Therefore, different numbers of pair excitations are prac- 
tically independent from each other |43l k4 l45l |46( . 

After solving the one electron problem (Eqs. |II3j ) 
and obtaining the single-particle eigenstates 4>i and their 
energies £j, the Hamiltonian for the interacting electrons 
can be written in second quantization as 

H = E Eicfci + i E VijkicfcfckCi. (5) 

i ijkl 



Here V^m's are the Coulomb matrix elements, 

v m = //«W«W 55 ^ r5r 

fe (r 2 )0i(ri)dridr 2 , (6) 

e(ri, r 2 ) is the dielectric function [3(|. We replace it with 
the dielectric constant e throughout the calculation. The 
method for computation of these elements is given in the 
appendix. The Hamiltonian for many-exciton complex 
can be written as [32l l4l| , 

i i ijkl 

+ J2 Xnkihfcfckhi + \Y1 V ifkl c t c t CkCl 

ijkl ijkl 

+ lj2 v w h t h t h * h t- ( ? ) 

ijkl 

The electron-hole exchange interaction elements, Xijki 
are defined by 

x m = ! [ ^(r^ (r 2 )— — e — r 

J J J 47re(n,r 2 ) • |n - r 2 | 

</> fe (ri)(^(r 2 )dridr 2 . (8) 

IV. EXCITON RECOMBINATION 

In order to calculate the photoluminescence spectrum, 
one needs to calculate eigenstates of both N exciton and 
N — 1 exciton systems. At low temperature, only the 
ground state and a few excited states of the N exciton 
system are required. However, in order to obtain the 
spectrum over a broad energy range, a larger number of 
eigenstates of the N — 1 exciton system has to be calcu- 
lated. In general, about 1000 — 2000 eigenstates of the 
N — 1 exciton system are required to cover transitions 
occurring in the s and p shells. 

Let us begin with recombination of non-interacting 
electrons. The momentum matrix element between an 
electron state 4> e = J2 n V'nWn and a hole state 4>h — 
E» ^"n is gi ven by 

ran m 

Here, ^v's are the envelop functions and the basis func- 
tions UnS are chosen as eight uncoupled spin-orbitals, 
i.e., \s T), \x T), \y ]),\z t), \s I), \x I), \y |), and \z [). 
If we neglect the contribution from the envelope-function 
part of the wave function J2 m (ipm\e ■ plV'm) SSL if can 
be further simplified as 

- «h£ T >-«IV£i>L 
(MPym - iP -M^) + ^ h vl \r si ) 

- (i>sMi)-{^i\r yl )], (io) 



where iP = (s\p x \x) — (s\p y \y) denotes the coupling 
between the conduction and valence bands. For circular 
polarization a + or er~, the momentum matrix element is 
then given by p± e = -^((cj) h \p x \<j> e ) T i(<j>h\p y \<f> e })- 

In the single-band effective-mass method, the Bloch 
functions for the heavy hole are Uh] — ^{\ x T) + Av T)) 
for j z = 3/2 and u H = ±(\x |) - i\y j)) for j z = -3/2. 
Hence, we have p^ e = (0 /lT |p~|0 eT ) = ( U/lT | -±= (p x + 

*Py)l s Diiphl^e) = -iPoi^hl^e)- It is straightforward 
to show that p~^ e = pt . 

The intensity of photoluminescence from the recombi- 
nation of one electron-hole pair in a iV-exciton state is 
defined as HHHl 

I a± (hv)= Y.^ E n)H^ C n-i\ P ^\ C n)\ 2 

f 

■SiEj, - E f N ^ - hv), (11) 

where C l N is the i-th eigenstate of the 7V-exciton system. 
Note that (C^ r _ 1 |P c ^ ± \C l N ) coherently sums all the possi- 
ble recombinations, therefore, the interference effect may 
play an important role. 

The probability function is defined as f{E l N ) = 
exp(-£^/rtT)/5]jexp(-£^/KT). The operator P~ ± 
describes all the possible electron-hole recombination, 
namely, 

P(T± = ^ ' Pimi^" c "ii (12) 
nm 

In the absence of magnetic field, we have I a+ (E) = 

lajE). 

V. RESULTS AND DISCUSSION 

We now illustrate our method by a calculation 
for a model structure of Ino.5Gao.5As/GaAs disklike 
SAQD characteristic of SAQD's grown using In- flush 
methodji?]. The dot has diameter 25.4 nm along the 
base and 2.3 nm height along the growth direction. The 
composition and dimensions of the dot are chosen such 
that its emission spectrum is similar to the one observed 
in the experiment of Raymond et al |49j . The strain dis- 
tribution is calculated by the continuum elasticity theory 
[Til IHof on a large cubic finite-difference mesh that has 
120 nm along each dimension and Dirichlet boundary 
condition on each side in order to ensure that the strain 
is fully relaxed. 

In Fig. |21 we show the probability density of the first 
six electron (in the upper row) and six hole states (lower 
row) calculated by the multiband k ■ p method. The cor- 
responding energy levels are plotted in Fig. [21 The ma- 
terial parameters used in the calculation are taken from 
Ref. Q . 

The circular symmetry of the single-particle states is 
found basically preserved due to the small shear strain 
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FIG. 2: Two-dimensional plot of the density of first six elec- 
tron (first row) and hole (second row) wave functions calcu- 
lated by the multiband k ■ p method. 



and weak piezoelectric potential in this intermixed quan- 
tum dot. Hence, the states in the conduction band and 
valence bands are seen to group into three shell, respec- 
tively. In the second shell, the shear strain induces a 
small splitting of 1.3 meV between the two p-like valence- 
band states. In the third shell, the splitting is about 
3 meV and the disklike geometry is responsible for the 
splitting between the two 3d states and the 2s state. It 
should also be noted that the ordering of the electron 
states in this shell is different from that of the hole states. 

In the single-band effective mass calculation, the four 
effective mass parameters are chosen to fit the spectrum 
to that obtained by the k ■ p method, giving ( in unit 
of free electron mass mo ) mjj = 0.060, = 0.070, 

m| = 0.27, and m\ = 0.30. It is noted that the electron 
effective mass in the dot is larger than the bulk value 
0.045 for Ino.5Gao.5As and approaches the value in bulk 
GaAs. The anisotropy in the effective-mass tensor in the 
valence bands, 0.30/0.27, is greatly reduced comparing 
with the bulk value, 0.29/0.074, i.e., the holes are much 
heavier in the plane perpendicular to the growth direc- 
tion. Similar findings that the effective mass of electrons 
in quantum dots exceeds the value in the corresponding 
bulk dot material and approach that in the bulk matrix 
material and the in-plane component of the effective mass 
of holes becomes much lighter have been reported [EOSl- 
In valence bands, the heavy hole and light hole are 
decoupled by the biaxial strain. For dots of small height, 
the biaxial strain is almost constant inside the structure, 
hence, the low-lying states in valence bands are mostly 
heavy-hole states. This is the reason why these states 
can be fitted by using single-band approximation. For 
thick dots, the band edges of heavy hole and light hole 
may cross each other due to the fact that the biaxial 
strain changes its sign inside the structure, which results 
in more light-hole components in the hole states in these 
dots. 

The energy levels calculated by EBOM differs from 
those by the k ■ p method, especially for the high-lying 
states. The shell separations by EBOM are smaller than 
those by the k ■ p method. For example, in the conduc- 
tion band, E s _ p (separation between s and p shell) is 24.8 
meV by EBOM and 27.3 meV by the k-p method, respec- 
tively. The averaged separations E p ^d between different 



5 




EBOM 



FIG. 3: Calculated energy levels of the Ino.5Gao.5As/GaAs 
self-assembled quantum dot by the eight-band k ■ p (KP) 
method and EBOM. Also shown are the energy levels fitted 
for the effective- mass approximation (EMA). 



duction band consist of less than 1% component from the 
split-off band while for valence band states it is less than 
5%. Hence, the mixture between spin up and spin down 
components in any of these states should be very small, 
i.e., they are polarized. 

However, in the absence of magnetic fields, all the 
single-particle states calculated by the k ■ p method and 
EBOM are doubly degenerate due to the time-reversal 
symmetry [lfj. Instead of having two degenerate states, 
the numerical calculation can only give one state from 
a random linear combination of the two polarized and 
degenerate states. Because of this degeneracy, most of 
the calculated single-particle states are found not polar- 
ized. By applying a small magnetic field (1 mT) along 
the growth direction, this time-reversal symmetry can be 
removed and polarized single-particle states are recov- 
ered. The eight-band k ■ p Hamiltonian is modified [5^ 
to include the effects of magnetic fields. For EBOM, we 
introduce Peierls phase factor [5^ to include the mag- 
netic field in the Hamiltonian. 



models are even larger, 30.6 meV from EBOM and 36.4 
meV from the k ■ p method. 

As shown in Fig. ^ the band structure predicted by 
the eight-band k ■ p method and EBOM match only in 
the region close to the T point. Although the lateral 
dimensions of the quantum-dot structure are large, it is 
very thin (2.3 nm) along the growth direction. Hence, 
the confined electron and hole states include components 
with large k z , which results in different energies for these 
states. For dots with larger height, there is very little 
discrepancy found between the two methods p4| . 

A. Polarization of single-particle states 

In the framework of the envelope function formalism, 
it is possible to separate the components of an envelope 
function into two groups when the spin-orbit interaction 
is not very strong. One group consists of components 
for spin up basis functions, \s f), \x |), \y f), and \z f), 
and the other consists of components for spin down basis 
functions. We define polarization of state <j> as 

p = J |^ T (r)| 2 + |^ T (r)| 2 + |^ T (r)| 2 + |^TW| 2 dr. 

(13) 

A state is polarized if cither p ss 1 (a 'spin' up state) 
or p k (a 'spin' down state). Apparently, there is 
little overlap between the polarized states with different 
polarization. 

A careful examination of the calculated single-particle 
states shows that all the single-particle states in the con- 



B. Electron- hole exchange interaction 

The ground state of a single exciton is a dark doublet 
separated from a bright doublet by the exchange energy. 
A dark (bright) exciton state is dominated by a configu- 
ration of an electron and a hole in their respective ground 
state with the opposite (same) spin. The bright doublet 
has a higher energy due to the electron-hole exchange 
energy. Because of the relatively large size of SAQDs, 
the electron-hole exchange interaction causes a very small 
correction to exciton states. It can be measured from the 
fine structure of single exciton recombination spectrum 
0. 

An accurate calculation of the electron-hole exchange 
will require knowledge of both electron and hole states 
and the dielectric function at a microscopic level |36l.l3^|. 
However, for SAQDs, the exchange energy can be esti- 
mated by using the multiband k ■ p theory or EBOM, 
where the electron-hole exchange interaction arises from 
the mixing between the conduction and valence bands. 
The calculation for our structure shows that the sepa- 
ration between the dark and bright doublets of a single 
exciton is 74.6 /ieV, which includes an electron-hole ex- 
change energy 63.2 fj,eV and the correlation effect. It 
gives a fairly good agreement with the value derived from 
the experiment 54] on similar samples, considering the 
approximation made in the theory, and the uncertainty 
in the dot size, shape and composition in the experiment. 
Because of the small contribution of the electron-hole ex- 
change interaction, we neglect it in the subsequent cal- 
culation. 
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C. Addition energies and hidden symmetry in 
multi-exciton complexes 

Due to the presence of quasi-shell structure in the 
single-particle energy spectrum, we chose the first 12 elec- 
tron and 12 hole states (with quasispins), which form the 
first three shells in conduction and valence bands, respec- 
tively, to build the multi-exciton configurations. 

In order to reduce the total number of configurations 
which grows factorially with the size of single-particle ba- 
sis set, we impose an additional constraint on the exciton 
configurations, i.e., the sum of the electron quasi-spins 
should be equal to that of the hole quasispins 0. We 
also apply a truncation according to the Hartree-Fock en- 
ergies of configurations in order to limit the total number 
of the configurations to less than 50, 000. 

Hidden symmetry [H, |U |3 03 is a good approxi- 
mation in a multi-exciton system with degenerate single- 
particle states and symmetric electron-electron, electron- 
hole, and hole-hole interactions. It predicts that the 
chemical potential i.e. the energy required to add an 
electron-hole pair to the system, is independent of the 
number of excitons. 

The symmetry between electron-electron, electron- 
hole, and hole-hole interactions is broken because the 
hole states are usually more confined than the electrons 
in the conduction bands. Hence, the hole-hole inter- 
action is generally stronger than the electron-electron 
interaction. The calculation by the multiband k ■ 
p (single-band effective-mass) method shows Vff n = 
14.4(14.9) meV, Vft^ = 16.7(19.4) meV, and V# e n = 
15.4(16.6) meV. The EBOM gives similar values, which 
are 13.9 meV, 16.3 meV, and 14.9 meV, respectively. In 
the case of either k ■ p or EBOM, the hole-hole interac- 
tion is stronger than the electron-electron interaction by 
about 15%. The hole states from the single-band calcu- 
lation are found more confined than those from the two 
multiband calculations. 

Fig. 0] plots the calculated addition energies for dif- 
ferent number of excitons. Both the k ■ p method and 
EBOM give similar result except that the EBOM predicts 
lower values. The result of the single-band effective-mass 
approximation is found very close to that by the k ■ p 
method. A clear plateau structure can be found associ- 
ated with the shell structure, which is an apparent signa- 
ture of hidden symmetry. The fluctuation in the addition 
energies within the same shell is not larger than that in 
the single-particle energies or the difference among the 
electron (hole)-electron (hole) interactions. The k ■ p 
method gives the largest fluctuation 1.4 meV in the p 
shell and 4.5 meV in the d shell, while EBOM gives 1.7 
meV and 2.7 meV, respectively. It is therefore seen that 
the hidden symmetry in our structure is well preserved, 
and not sensitive either to the splitting of p and d shells 
or the asymmetric interactions. 
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FIG. 4: Addition energy spectra /n(7V) = E g (N) -E g (N- 1), 
calculated by the single-band effective-mass approximation 
(triangular dots), the multiband k -p method (circular dots) 
and EBOM (square dots) for the Ino.sGao.s As/GaAs disklike 
quantum dot. 
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FIG. 5: Excitonic emission spectra calculated by the single- 
band effective-mass method. 



D. Emission spectra 

Figs. 03 El and [7] show the emission spectra calculated 
by the single-band effective-mass method, the eight-band 
k ■ p method, and EBOM, respectively, for up to six ex- 
citons. As the k ■ p method and EBOM are shown to 
give very similar results except for overall shifts in tran- 
sition energies, we will confine our attention to analyzing 
the difference between the results by the single-band and 
the multiband k -p calculations. It should be mentioned 
that the spectra in Figs. and El are plotted in different 
scale, only the relative intensity between the spectra in 
the same figure is relevant. 
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FIG. 6: Excitonic emission spectra calculated by the multi- 
band k ■ p method. 
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FIG. 7: Excitonic emission spectra calculated by EBOM. 



The difference between the single-band and multiband 
calculations concentrates in the s shell as the emission 
peaks in the p shell are found similar between Figs. 
and El Compared with the single-band calculation, the 
band-mixing effect enables more configurations from the 
multiband single-particle states to be coupled through 
the Coulomb interaction. It results in more possible yet 
weak transitions in the low energy end of the spectra, as 
are found in Figs. El and [7| when the number of excitons 
is larger than 4. 

This is illustrated in Fig. 00 In the emission spectra 



from the three-exciton complex (3X), all the three meth- 
ods give one strong peak in the p-shell region and exhibit 
different structure in the s-shell region. The single-band 
calculation gives one peak with high intensity and three 
other small peaks in the s-shell region while both the k-p 
method and EBOM show five peaks with visible intensi- 
ties. 

The initial state of the recombination from 3X is its 
ground state and the final states associated with peaks 
found in the s-shell region are excited biexcition states 
|34| . As there is little difference among the 3X ground 
states calculated by different methods, it is the excited 
biexciton states that account for the different structure 
in the emission spectra. 

In the single-band calculation, the emission peak with 
high intensity in the s-shell region is associated with three 
biexciton states where both the two electrons and the two 
holes are in a triplet configuration (one in the s shell and 
the other in the p shell). In these states, the total spin 
S(= I) and its z-component S z of the two electrons are 
the same as those of the two holes. As the spin-orbit 
interaction is absent in the single-band calculation, these 
three biexciton states of different S z are degenerate and 
give rise to only one peak in the s-shell region. 

In the presence of spin-orbit interaction which is taken 
into account in both the k-p method and EBOM, the de- 
generacy among these three biexciton states is partially 
lifted, i.e., the biexciton state of \S Z \ = I has a different 
energy from that of Sz = 0. It gives rise to two splitted 
peaks in the s-shell region, as illustrated in Fig. [SJ The 
three other smaller peaks are associated with those ex- 
cited biexciton states in singlet-singlet or singlet-triplet 
configurations, which are less affected by the spin-orbit 
interaction, and can be seen in the all spectra. 



VI. CONCLUSIONS 

In conclusion, we have presented a multiband mi- 
croscopic theory of many-exciton complexes in self- 
assembled quantum dots. Three methods: single-band 
effective-mass approximation, the multiband k-p method, 
and the tight-binding-likc EBOM, are used to obtain 
single-particle states. We expand the many-body wave 
functions of N electrons and N valence holes in the ba- 
sis of Slater determinants. The Coulomb matrix ele- 
ments are evaluated using statically screened interaction 
and the correlated 7V-exciton states are obtained by the 
configuration interaction method. We apply the theory 
to the excitonic emission spectrum in InAs/GaAs self- 
assembled quantum dots and successfully compare the 
results of the single-band effective-mass approximation 
with those obtained by using the of k-p and tight-binding 
methods. 
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FIG. 8: Emission spectra from three-exciton complex calcu- 
lated by different methods. 
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APPENDIX A: COULOMB MATRIX ELEMENTS 



these localized atomic orbitals and replace ri — r2 with 
Ri — R2 in Eq. [fjl we have 

R.!,R.2 m 

^^£^R 2 )^R 2 ).(A2) 

1 1 n 

We further transform a three-dimensional (3D) wave 
function ip(R) into a column vector tp(r) by mapping 
the 3D variable R onto a one-dimensional index r, a 
six-dimensional integral can be converted into a vector- 
matrix multiplication, 

V ijkl = £ <8 Vwf -U-(52 ^*jn ® i>hn), (A3) 

m n 

where (£> is the direct multiplication (element by element) 
operator between two vectors. U is the matrix with el- 
ements U{r\,r2) = 4 7re .|R 1 _R 2 | ■ The diagonal elements 
can be obtained by the integration of 1/R over a unit 
cell. 

In order to use the optimized BLAS (Basic Linear Al- 
gebra Subprograms) library [56j, the above formulation 
can be further vectorized as 

J = $ T •[/•$, (A4) 

where & r ,{ij} = ^'iraMV'jnM is a matrix contain- 
ing all the possible pairs of two-particle wave functions. 
Due to the large dimension of matrix U , we make use 
of domain decomposition in the numerical calculation to 
divide it into a number of smaller matrices and sum up 
the result of all the individual multiplications. 



In this appendix, we describe an efficient approach to 
calculate Coulomb matrix elements numerically. In the 
configuration interaction method, the properties of the 
system are given by the single particle spectrum and by 
the Coulomb matrix elements defined as two-electron in- 
tegrals (see Eq. |H]). One possible way to calculate these 
six-dimensional integrals is to first solve Poisson equa- 
tion, than calculate reduced three-dimensional integrals 
|55j . The calculation is repeated for each integrals. For 
the calculation involving 12 electron states and 12 hole 
states, the total number of integrals is almost ten thou- 
sands. Here we propose an algorithm that does not re- 
quire calculation for each integrals. 

Within the envelope function formalism, the wave 
functions of single-particle states are expressed as a linear 
combination of Bloch sums, 

<Mr) = -^5353^„(RK,(r-R), (Al) 

* n R 

where ipi n is the n-th component of the envelop func- 
tion fy, u n (r — R) is the n-th atomic orbital localized at 

unit cell R and l/V^V"X)n u n( r ~ R-) is the correspond- 
ing Bloch function. If we ignore the contribution from 



APPENDIX B: COMPANION CONFIGURATION 
AND ADDITIVITY RULE 

In this appendix, we point out how to use the addi- 
tivity rule to construct configurations for multiexcitons. 
As a multi-exciton complex contains two different par- 
ticles, electrons and holes, the total number of possible 
configurations is much larger than for electrons or holes 
separately. To circumvent this difficulty, we use the fol- 
lowing algorithm for construction of multi-exciton con- 
figurations. 

The many-body Hamiltonian matrix constructed from 
the CI method is a sparse matrix. For a given config- 
uration, there is only a small number of configurations 
interacting with it, which are named as its companion 
configurations. Let us first define the distance between 
configuration C\ and Cj, ||C,,Cj||, as the total number 
of single-particle states that the two configurations differ 
by. It is apparent that 

(Ci\H\Cj)=0, if ||Ci,C,|| > 2. (Bl) 

An exciton configuration consists of a part for elec- 
tron^) and the other part for hole(s), i.e., Cf x = 
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{Cf,C^}. The distance between two exciton configu- 
rations, Cf x and Cj x , can be easily calculated by the 
additivity rule, namely, 

d% = \\cr,q x \\ = \\cf,q\\ + \\c^c^\\ 

IK, • IK, (B2) 

One can calculate the distance matrix D e and Dh for 
the electron and hole configurations, respectively, and 
then obtain the matrix D ex by using the additivity rule. 



When the number of single particle states (cither elec- 
trons or holes) is large, a cutoff is necessary to be applied 
to the total number of electron or hole configurations. 
Depending on the memory available for the computation, 
it is set to be between 5000 and 10000. Once the distance 
matrix D ex for the exciton configurations is calculated, 
it is straightforward to construct the configuration inter- 
action matrix as the positions of all the none-zero matrix 
elements are known. 
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